Light Propagation and Paired Superradiance in Coherent Medium 
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ABSTRACT 



The problem of light propagation of frequency corresponding to half of the energy difference 
O ■ between a metastable excited state and the ground state of atoms is examined, and solved for 



coherent medium by analytic means. We demonstrate that the non-linear system of Maxwell-Bloch 
equation for the effective model of the A— type three levels is integrable in the mathematical sense. 
Analytic solutions thus obtained describe pulse splitting accompanied by compression, indicating 
a kind of non-linear instability of propagating pulses. The instability is eventually terminated by 
coherent two photon emission (called paired superradiance or PSR in short). These results are 
displayed by numerical outputs for visual understanding, as well. It is further shown that the 
integrable system allows a new class of soliton solutions. Solitons, implying the phenomenon of seff- 
induced transparancy at non-resonant frequencies, are stable against PSR. One of our goals of the 
present work is construction of a calculable theoretical framework for PSR rates associated with a 
trigger pulse propagation, which is achieved by combining analytic results with perturbative methods. 
PSR photon spectrum and its rate oc(target number density) 2 , along with their time structure, are 
clarified this way. These results may open a new path for interesting technological applications such 
as quantum entanglement and for solving the remaining problems of the still mysterious neutrino. 
Some basic strategy for realistic experiments of PSR detection and soliton production is also outlined. 
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I Introduction 

Cooperative phenomenon known as superradiance (SR) pQ may have more dramatic effects when 
applied to the forbidden transition: the decay rate via two photon emission may be enhanced to 
a macroscopic level without the wavelength restriction [2]. This is in contrast to the single photon 
SR enhanced decay, with the coherence region limited by the wavelength 2 . We hereafter call the 
macro- coherent two photon emission [2] as paired superradiance (PSR). An example of candidate 
atoms for PSR is the first excited D-state of two electron system of Ba (its levels shown in FigflT])). 

In the course of establishing a firm theoretical formulation for PSR, we came accross a basic 
reference [5] (and quite possibly many related papers unknown to the present author) in which the 
Maxwell-Bloch or rather a Maxwell-Schrodinger equation for an effective two level problem is derived 
and some numerical simulation based on this equation is performed. 

The purpose of the present paper is four-fold; (1) to add to the literature of this field, fully 
integrated analytic solutions of the relevant non-linear system, which describe pulse splitting and 
compression, (2) to present a new class of soliton solutions of quantized pulse area (the terminology 
to be defined) that describe the phenomenon of self- induced transparancy(SIT), (3) to give the 
PSR rate and its spectrum, both time dependent, along with its relation to soliton formation, (4) 
to provide the signal to the noise ratio (S/N) for radiative neutrino pair emission RNPE (a new, 
proposed experimental method for the precision neutrino mass spectroscopy, [I], [5]). 

Our analytic solutions describe the multiple pulse splitting accompanied by compression of prop- 
agating pulses in a long target medium, when the input pulse is strong enough. The pulse strength 
is made quantitative by our own concept of the pulse area, a product of integrated field flux and 
the coupling strength to the medium. The phenomenon is a highly non-linear coherent effect of the 
whole system of target atoms and fields. 

The multiple pulse splitting in a long target may further be interpreted as a process towards 
formation of many solitons, if the effect of PSR is included. Indivisual split pulses in medium 
become increasingly sharper prior to PSR. The ever sharpening pulse has an ever increasing energy 
density of fields and is unstable against the physical process of PSR. It must become stabilized by 
PSR emission, eventually resulting in formation of solitons, objects stable against PSR. This intuitive 
picture has been supported by extensive numerical computations performed by the present author, 
some of which are shown below. We would like to convince even the uneducated reader of this simple 
picture of what occurs in an ideal envirornment, the infinitely long coherent medium. 

Solitons of the quantized pulse area and of arbitrary velocities have much simpler analytic forms 
than in the case of two level system [5J. We give explicit formulas and exhibit their pulse shape in 
a number of figures. Moreover, SIT occurs at non-resonant frequencies, hence might be more useful 
in technological applications such as communication by light and quantum entanglement: excellent 
candidates are (J=0 — > J=0 ) PSR transitions in alkhali earth atoms. 
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Our general method also allows to discuss residual interaction between two solitons, which turns 
out attractive, suggesting existence of bound states of two solitons, presumably even a possibility 
of field condensate. Implication of the attractive nature of force shall be addressed in forthcoming 
work. 

We next compute, based on perturbative methods, PSR rates in which one of the photons is the 
forward going pulse component. In PSR, two photons are, almost exactly, back to back emitted and 
have the same energy, or only this configuration in the two photon phase space is macroscopically 
enhanced. Calculated rates are more than adequate for detection, and even give a hope of PSR 
measurement in ion traps where one can expect to store a total target number of ions only as large 
as 10 6 . 

The time structure of our triggered PSR is complicated. PSR is expected to occur rapidly and 
violently, most dramatically just prior to the soliton formation (and at the time of its artificial de- 
struction). After soliton formation PSR stops, and propagating pulses become stabilized as solitons. 
An ideal, and our own favorite, method of observing PSR is creation of as many as possible solitons 
and their subsequent controlled destruction, giving the largest PSR rates at the instants of creation 
and destruction. We would also like to stress, for PSR detection and soliton production, great ad- 
vantages of non-resonant frequency of the trigger laser, which is non-destructive to target atoms and 
makes experiments easier. 

In any of these phenomena, irradiation of the triggering laser at the doubled wavelength is of 
vital importance to us, and we must fully understand the problem of light propagation in medium, 
which is done below. In the literature [I] that deals with SR of the two level atom, SR initiated 
by quantum fluctuation has extensively been discussed. The initial setting that interests us most 
in the present work is considerably different from this initial condition. PSR initiated by quantum 
fluctuation is weaker and moreover it is our main intention to utilize benefits of the triggering laser. 
Thus, our problem is more akin to the triggered or induced SR in the two level problem [I] . Related 
to this, the time delay observed in fluctuation initiated SR is absent in our triggered PSR. Despite 
of all these, we shall briefly mention interesting features of almost trigger-less PSR under a weaker 
trigger; excellent experimental signatures along with a rough estimate of its rate. 

In most of our analysis below, we ignore relaxation processes, due to that the enhanced PSR is 
very fast: this seems a legitimate approximation under a wide range of circumstances. Moreover, the 
coherence we need for PSR is not a stationary state of target atoms, but rather it is a dynamically 
developed (with time) state of the whole system of target atoms plus fields due to the non-linear 
interaction between the two system. In actual experiments, the coherence is dynamically generated 
during a short time interval via a series of physical processes of excitation and trigger by lasers. If the 
time of coherence development is shorter than a typical relaxation time of relevant target state, one 
could achieve a well prepared state for PSR. In this sense what is needed for a positive measurement 
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is, in addition to the fast PSR rate, a fast preparation of coherence development, whose realization 
is left to a challenge for experimentalists. 

At a fundamental level the soliton formation may give rise to a controlled measurement of RNPE, 
0]) 0, because solitons are stable against two photon emission, regarded as a crucial background to 
RNPE. Solitons are, on the other hand, unstable against RNPE, hence soliton formation enhances 
the signal to the background ratio in favor of RNPE. If the enhanced RNPE rate is larger than 
(spontaneous decay) rates of the next leading order QED processes (usually much smaller than 1 
msec -1 order), then the PSR background suppressed RNPE (due to two effects discussed in the 
present work; mismatched trigger frequency to PSR and the soliton formation) is measurable by well 
controlled experiments. 

The present paper is organized as a collection of (hopefully) compact sections containing many 
numerical outputs, and two long mathematical appendices: section headings are II Effective two level 
model, III Maxwell-Bloch equation and equation for the tipping angle of Bloch vector, IV Construc- 
tion of analytic solutions, V Soliton solutions, VI Theory of PSR and its relation to soliton formation, 
VII Outlook for RNPE, VIII Appendix I Derivation of effective two level model, IX Appendix II De- 
tails towards construction of analytic solutions. We present numerical outputs, presumably more 
than necessary to the educated reader, to help even the uneducated reader easily understand physics 
behind analytic and numerical results. Appendices give some lengthy details of derivation omitted 
in the main body of the text. 

It is hoped that this work helps experimentalists to design clever methods of detecting and 
measuring PSR and of creating optical soliton of our type. 

Throughout this work we use the natural unit so that h = 1 ,c = 1. 



6s6p 'p, 




6s 1 'S„ 



Figure 1: A— type low lying levels of neutral Ba atom. An excitation scheme to 1 L*2 is also illustrated. 
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II Effective two level model 

In this section we follow [3] to sketch derivation of the coupled non-linear partial differential 
equations of Maxwell-Bloch (MB) type, with a slight modification of notation and correction of 
mistakes at detuned frequencies. The essence is to derive an effective two level model for atoms of 
A— type levels such as Fig(jTj). 

The spirit of this model is a long time average, or a truncation of past memory effects, in the sense 
of Markovian approximation, and slowly varying envelope approximation (SVEA). These should be 
an excellent approximation in our light propagation problem in medium, because this light wave 
at the doubled wavelength is non-resonant to target atoms, its energy being too far away from the 
frequency to excite to the first atomic level above the ground. 

The model is summarized by the Schrodinger equation for two amplitudes of two lower levels; a 
S-state, 1 5'o (its amplitude given by c g ) and a D-state, 1 D 2 (given by c e ) in the case of Ba atom. It 
has the structure of two x two quantum system of an effective Hamiltonian H; 




e -i{^-Ee 9 )t^ E *Y 



(1) 



(2) 



Here E (x, t) is the slowly varying envelope of fast oscillating field (to be multiplied by e luJt ). We only 
discuss one dimensional propagation problem, assuming the axial symmetry around the propagation 
axis, taken here 

Coupling strength here is given in terms of electric (or magnetic) dipole moments between the 
upper level P-states in the Ba case (the leading candidate being 1 Pi), denoted by dj e , dj g ; 

j2 p j2 rp 

_ 2 V a je^je -oV j9 j9 (3) 

j je M j ^jg M 

j j 
E c = ±{E jg + E je ), Su = u-^. (5) 

El or Ml moments dji are related to experimentally measured (or theoretically calculated in some 
unfortunate cases) decay rates 7^ by d 2 ^ = 3ir r yji/Ej i with Eji the known energy difference between 
two levels. Higher levels denoted by \j) are assumed connected to two lowest levels, \e) and \g), by 
strong transition elements. It would often be sufficient to consider a single level for \j) of strongest 
coupling. As explained in Appendix I, one may assume to a good approximation at u> ~ E eg /2 the 
symmetry fi eg = fi ge , which we take in the rest of the text. 

Matrix elements in H, eq.([2]), of our effective A— model are due to the Stark effect. It is important 
to recognize that this shift contains two powers of propagating field Eq, its absolute value being the 
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power or the flux of propagating field. The magnitude of Stark shift is not negligible for strong laser 
irradiation: two diagonal elements are of order, 6GHz for the D-state and 16 GHz for the S-state, 
and 2.1 GHz for the magnitude of off-diagonal elements, assuming the laser power of 10 6 W mm -2 . 
Derivation and neglected terms of this Hamiltonian structure are discussed in our Appendix I. 



Ill Maxwell-Bloch equation and equation for the tipping angle of Bloch vector 

The next step is to derive a macroscopic coupled set of equations for polarization of medium and 
propagating field. Throughout this work we take the continuum limit of distributed target atoms, 
and assume an axial symmetry around the axis of pulse propagation, which reduces our problem to 
1 + 1 (one time t and one space coordinate x) dimensional field theory. Polarization of medium is 
defined by a Bloch vector of 3 components (a part of density matrix elements R a p in the notation of 
Appendix I); 

R 1 = m(c*c e e l " - c* e c g e^) , (6) 
R 2 = -n(c* g c e e 111 + c\c g e~^) , (7) 
R 3 = n(\c e \ 2 - \c g \ 2 ) , (8) 
r] = (2u- E eg )t - 2kx + 2(p . (9) 

Here n(x) is the number density of targets assumed at rest, and we took for simplicity a constant 
density n(x) = n within a medium of finite length, < x < L. The function ut — kx + ip(x, t) is the 
main phase part of the propagating field e'^Eo. 

From the Schrodinger equation we derive the Bloch equation for polarization components, 

^Ri = ^^\E 2 \R 2 + f^\E 2 \R 3 , (10) 



dt 1 4 

9 R 2 = -^^\E 2 \R 1 , (11) 



dt 

!* s = -f|£Sl«, (12) 
From this set of equations one has a conservation of the magnitude of the Bloch vector, 

|( J R 2 + J R 2 + J R 2 ) = 0, (13) 

hence R 2 (x,t) is time independent, which can be taken the squared number density n 2 (x). 
A linear combination of R 2 and R%, 

= R 2 -jR« (14) 
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is conserved in our system of differential equations, and one may set this vanishing without any loss 
of generality, giving the condition R 2 = 7.R3. The Bloch equation is thus effectively reduced to 

!^i = ^(l + 7 2 )|£o 2 |^3, (15) 
§-R 3 = -^\EZ\ Rl , (16) 

with a constant of motion, Rf + (1 + 7 2 )-R 2 = n 2 (x). For the Ba D-state, I7I ~ 2.3. 

Strictly, the Bloch vector has 4 components, but to a good approximation the other component 
is time independent. We shall ignore this complication, relegating some explanation to Appendix I. 

To proceed further, it is convenient to introduce the angle function 9(x,t), called the tipping 
angle, and the constant of motion B = ±n(x) by 

R l = Bsin9, R 3 = = cos9. (17) 

V 1 + 7 2 

The Bloch equation is then equivalent to a relation between 9 and the field power \Eq\\ 

d t 9 = u R , u R = fx\E 2 \, (18) 

J I + 7 2 1 / 

P-= 2 fJ,ge= 4\^ 99 ~^ e ) 2 + 4/i2 e , (19) 

with /1\Eq \ ~ 3G¥Lz\Eq |/(10 6 Wmm~ 2 ) for the Ba D-state. The tipping angle 9 is thus an integrated 
flux of \Eq\, with a weight given by the strength of Stark shift ft, 

/t—x 
dy\E%(x,y)\. (20) 
-00 

9{x, 00) is called the pulse area, while 9(x,t) the area function in the present work. The relation 
( 120 p . containing the squared amplitude \Eq\, is different from the corresponding one in the two 
level problem [TJ in which the field E Q itself appears in the integrand. As an illustration, the area 
9(x = 3cm , t = 00) is plotted as a function of input laser power in Fig(j2J), which we need to find out 
necessary power values for subsequent computations on the Ba D-state. Moreover, the pulse area of 
2tc is an important unit for quantized solitons, a subject fully discussed in the following sections. 

It is important to distinguish two cases of different signs of the constant B. If B = n > 0, the 
population difference R3 is positive for small 9, which means that there are more atoms in the excited 
state than in the ground state. If B = —n < 0, there are more atoms in the ground state than in 
the excited state. From obvious reasons, we use the terminology similar to the case of the two level 
problem: the terminology of amplifier is used for B > and the absorber for B < 0. Needless 
to say, this distinction is interchanged by a redefining transformation 9 — > 9 + ir. Nevertheless, 
it is customary to take the limit value of the tipping angle, 9(— 00) — > 0, hence we stick to this 
terminology. 
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Input power vs area for Ba 

area/2 Pi 




W mm A -2 



2xl0 6 4xl0 6 6xl0 6 8xl0 6 



Figure 2: Pulse area of Ba D-state vs the input pulse power. Quantized areas of an integer number 
x27r can give rise to this number of solitons. Pulse duration 3 ns assumed. 

Polarization Ri of medium is related to the field Eq by the Maxwell equation. We discuss right- 
moving field given by \E (x, t)\e~ lul<yt ~ x \ Under SVEA, the evolution of envelope \E (x,t)\ is given 
by 

(«9i + <9 x )|£ | 2 = w/v|£ol#i- (21) 

The coupled system of non-linear partial differential equations, eq. (fT5|) . eq. (fT6]) and eq. (l2Tj) . forms 
the fundamental equation of field propagation in medium, which may be called the Maxwell-Bloch 
equation for our effective A— model. 

The Maxwell-Bloch equation becomes an equation for the single variable 8(x,t), using eq.( TT71) ; 

(d t + d x )d t 9 = ±a sin 9d t 8 , a = ujfi ge n . (22) 

Our basic problem is then to solve this single, non-linear equation under an arbitrary initial and 
boundary data. Difference from the two level problem is in the RHS term; in the two level case there 
is no dtO and only term cx sin 6, which is a potential term in the ordinary sense. This is the familiar 
sine-Gordon equation for propagation of pulses and SR at the resonant frequency [TJ. 

A simplest form of relaxation may be introduced into the envelope equation fl2TT) as a friction of 
the form, oc ftl-E^I- This leads to a modified equation for 8; 

(d t + d x )d t 6 + Kd t 6 = ±a sin 0d t . (23) 

Both forms of equation, eq. (1221) and eq. (1231) . may be integrated once, introducing an arbitrary function 
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A(x) of space coordinate; 

(d t + d x )9 + k9 ± a cos 9 = A(x) . (24) 

There are more complicated forms of relaxation such as the inhomogeneous Doppler broadening 
in gas [TJ, which has to be treated separately. 

When the initial data obeys the condition, |A(a;)| > a, the envelope may grow unlimitedly. This 
is a situation we do not discuss as our physics problem. With \A(x) \ < a, we set A(x) = ±a cos 6i{x). 
For the dissipationless case of k = 0, the once-integrated equation becomes 

(dt + d x )9 ± a(cos 9 - cos 9,) = . (25) 



IV Construction of analytic solutions 

It would be instructive first to discuss homogeneous solutions without spatial x dependence. 
The ordinary differential equation with the simplest relaxation included is given by 

cP6 . „. d6 

- + ( KTas i„ e) - = 0. (26) 

This describes a dynamical system of a fictitious pendulum (its angle location given by 9) with the 
friction term varying periodically with 9, k =f a sin 9. There is no force acting in the ordinary sense. 
Equation (|26l) may be once integrated to give the velocity d9/dt = —(k,9 ± acos#), after a suitable 
redefinition of 9 variable. The function of 9, —(k9/ol + cos9), is plotted in Fig(|3]) for large and 
small frictions, to illustrate importance of the magnitude k/oc. The region of d6/dt > alone is 
allowed by the energy positivity \Eq \ > 0, and this restricts the 9 region to be given by an inequality, 
— (k9 + a cos#) > 0, forming islands of 9 values nearly periodically separated for small n/a. For the 
dissipationless case of k = 0, there are infinitely many pieces of such islands, while these regions 
are limited in number for a finite k. This dynamically allowed number of islands gives the possible 
maximal number of produced solitons, as clarified later. 

The magnitude of a for Ba is ~ 0.72 GHzn/(10 16 cm -3 ). Smaller values of friction, e.g., less 
than 10~ 2 a (the factor 10~ 2 somewhat arbitrarily chosen), imply relaxation times larger than ~ 
140 ns (n/10 16 cm -3 ) -1 . We may also compare these values with the Doppler broadening in Ba 
gas, Ajj ~ 680MHz(T/300i^) 1 / 2 (with T the temperature). The requirement < a gives T < 
340i^(n/10 16 cm -3 ) 2 . We assume below that k, Ac and all other relaxation rates (including those in 
solids) are much smaller than a, and consider the dissipasionless case. 
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small(O.OOl) and large(O.Ol) relaxation 




Figure 3: Allowed islands of the tipping angle given by positive ordinate values. Two choices of 
friction, k/o, = 10~ 2 (in blue) and 10~ 3 (in purple), are compared. Smaller frictions give more 
islands, giving a formation chance of more solitons. 



The explicit solution without space coordinate dependence is, in the dissipationless case of k = 0, 

cos 6»o cosh (a(t- t ) sin 0„) - 1 

cos 9 = — + \ • a \ T~ ' ( 27 > 

cosh (a(t — to) sm 9 ) — cos 9q 

Ur _<^__ a sin2 / 2 gN 

dt cosh (a(t — t ) sin 9 ) — cos 9 ' 
^ _ ^ cos g cosh (a(t - t ) sin 6> ) - 1 ^ 2g ^ 
cosh (a (t — t ) sin ^o) — cos 0o 

This solution describes a dynamic motion of 9 that starts from 9$ at t = — oo, reaches 7T at t — to, 
and ends finally at 2tt — 9 . 

The method of deriving more general solutions of the full partial differential equation is to let 
parameters to , $0 here to depend on other variables of t — x , x. Details of this construction are given 
in Appendix II. 

Simplest solutions that describe the initial target state either fully in the ground state (denoted 
by ^ and called the absorber) or in the excited state (denoted by ( e > and called the amplifier) are 
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given by 



\E 2 Q (x,t)\ 



4(t-x) 



\E 2 (x,t)\ 



Ri 9) ~ -n 



1 - ax sin # + (ax) 2 sin 2 (0/2) 
cos 9 — ax sin 6 + (ax) 2 sin 2 (#/2) 
1 — ax sin 9 + (ax) 2 sin 2 (6/2) 

(g) _ 6g(t — X) 



1 + ax sin 6 + (ax) 2 sin 2 (0/2) ' 
cos 6 + ax sin 6 + (ax) 2 sin 2 (#/2) 
1 + ax sin + (ax) 2 sin 2 (#/2) 



6(t — x) = n 

e^(t — x) here is the input pulse intensity. 



t—x 



dyel(y) . 



Field power in amplifier: area= 1.04*2 Pi, l=3cm 



W mm A -2 




(30) 
(31) 
(32) 
(33) 
(34) 



Figure 4: Compressed pulse (in blue) compared with incident pulse (in purple). The time delay 
(measured in difference of central locations of the two peaks) common in SR is absent here. The 
absolute value of abscissa coordinate is meaningless. Number density 10 17 cm~ 3 , target length = 3 
cm, laser of power ~ 2tt pulse, duration 3 ns, and the initial angle 6 = n/4 are assumed. 



The most general form of solutions is given in Appendix II, eq. (1141j) ~ eq. (11421) . They are written 
in terms of an angle factor 6q and the same area function 6(t — x) as eq. f)34p . The angle 8 describes 
the initial state of targets. The limit 6q — > or it, describing the fully excited or the fully ground 
state, gives the solutions above. A finite value of 6q gives a kind of the dark state [8], a measure of 
quantum mixture of the metastable excited and the ground state. 
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Field power in amplifier: area=4.01*2 Pi, l=3cm 

W mm A -2 




Figure 5: Multiply split pulses in medium over the entire region of incident pulse. Number density 
10 17 cm -3 , target length = 3 cm, laser of power ~ 4 x 2tt pulse, duration 3 ns assumed. The initial 
angle #0 = tt/4. 



Two important parameters that describe solutions of eqs. ( |30l) ~ ( 1341) are a and the dimensionless 
constant K in the area function 6, and they are given by 

where 5 is the pulse duration. For the Ba D-state, they are numerically 

1 n 1.53a; . , 

a ~ 0.052cm- 1 — -, 37 

10i6 C m- 3 1.53 2 - (u - 0.707) 2 ' v ; 

ns 2.24 2 - 10 6 Wmm- 2 ' 1 ' 

with w to be given in the eV unit. 

For targets of short length, we may globally characterize the pulse modification in terms of the 
gain and the loss. Relative to the input pulse power (irradiated at x = 0) the transmitted pulse at 
the end of target (placed at x = L) has the gain or loss factor, 

^ ^— ■ (39) 

l±aLsinfl + (aL) 2 sin 2 fl/2 

We may think of sufficiently short pulse so that remaining polarization and field inside the target is 
effectively described by those values at an infinite time. The area function 9{y) is then replaced by 
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its value at time infinity 6(00) = 9^, hence the input pulse area 



Ooo = V / dyel(y) , 



(40) 



j —00 



may be used. The gain or the loss factor of the transmitted pulse is thus given by 



G 



1 



(41) 



1 ± ctL sin + (aL) 2 sin 2 9^/2 



A positive gain G > 1 requires, for < 9^ < 2tt, 



aL < =f2 cot — . 

2 



(42) 



For targets of longer length, the situation is more complicated: the pulse splitting occurs along 
with compression, as illustrated in Fig(jl|)~ Fig®. Both of these phenomena may be interpreted as 
signatures of instability of propagating pulses in medium. Time profiles of propagating pulses, along 
with the incident pulse without medium effects, are shown in Fig((l]) for the initial area of ~ 2ir and 
Fig© for the area ~ 8tt; both close to integer multiples of the quantized unit 2ti (the unit becoming 
important in discussion of solitons below). One does not observe delayed pulses, which might have 
been visible as shifts in two central peaks of pulses with and without medium. This is in contrast 
to the single photon SR initiated by quantum fluctuation, where the delay is universally present. 
If the target is long enough and the laser type is close to CW (continuous wave, and not pulse), 
there may be many split pulses within medium. In Fig(j6]) we show splitting into 8 well separated 
pulses. The number of separated pulses is ~ 6 l 0O /27r. Spatial profiles vary with time, as illustrated in 
Fig((71). The ever sharpening pulses are observed, as shown in Fig®, but in actual situations paired 
superradiance, not considered here, occurs and these sharpened pulses are expected to be rounded 
off by PSR and become objects close to solitons for which the pulse shape is unchanged with time. 

In summary, the pulse modification is described by two parameters; aL ,9^. The larger aL is, 
the larger modification occurs, while 9^ gives a measure of pulse splitting; the number of split pulses. 
In the Ba example, aL ~ 0.024nL/10 16 cm" 2 and 9^ ~ 4(5/ns)(eg/10 6 Wmm~ 2 ), both at u = E eg /2. 

We now discuss interaction of two pulses, assuming that both pulses propagate in the same 
direction. The calculation is possible because we solved the one-mode basic equation in terms of 
arbitrary initial data which can be a sum of two pulses of well separated envelopes at initial times, 
one of them catching up the other. The influence of 2nd pulse on 1st pulse is determined by using 
the area sum of two pulses, for the case of 9 = 0, 9i(t — x) + O^if — x), which replaces 9(t — x) 
in the formula, eq. (13"UlW (13"3l) . Suppose that two pulses are well separated by a distance X. Near 
the center of 1st pulse, we may take t — x ~ 0, then the effect of 2nd pulse is given by a pulse 
intensity modification factor A (to be multiplied by the energy density e 2 (0); the incident pulse 
intensity). This factor may be interpreted to define an effective potential V by using the relation 
A = 1 — V(X). This is like defining the potential V of relativistic particle from the Lagrangian L, 
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Pulse splitting: area= 8.76*2 Pi 

W mm A -2 




Figure 6: Splitting into many pulses in a long target. Laser (profile in purple) of intensity corre- 
sponding to 8.76 x27r area, time duration 50 ns and 5 GHz width is irradiated to a target of Ba 
number density 10 16 cm~ 3 and length 100 cm. 



Time varying profile: area=4.16*2 Pi, n= 10 A 16 cm A -3 



W mm A -2 



0.1 - 



0.001 




50 100 150 200 250 300 



Figure 7: Time evolution of field profile at 4 different times (blue — > purple — > brown — > green in the 
time sequence) spaced by 0.5 x 3 ns (pulse duration), 9 = ir/4. Input area close to 4 x 2n is taken. 
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Ever sharpening pulse: area=1.04*2 Pi, n= 10 A 16 cm A -3 

W mm A -2 



2.5 x10 s - 



2.0 x10 s - 



1.5 x10 s 



1.0x10 s - 



5.0 xlO 7 
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200 



250 
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cm 



Figure 8: Ever sharpening pulse, at 3, 6, 9 ns after the pulse input. 

using L = —mc 2 \Jl — j3 2 — V. For 6q = the potential becomes of the form (=p corresponding to 
amplifier/ absorber) , 

V(X) = 1 = -y- ^ , (43) 

1 =F axsin(0i(O) + 6 2 {X)) + (ax) 2 sin 2 (^(0) + 9 2 (X))/2\ 

where x is the distance 1st pulse propagated in medium (2nd pulse is at x — X in medium). This 
formula is of generic validity, and we shall apply it to the problem of two soliton interaction in the 
next section. 



V Soliton solutions 

Existence of stable solitons is anticipated from the topological reason when the pulse area takes 
quantized values of an integer times 2n. This situation is similar, but not identical, to the self- induced 
transparancy (SIT) in the two level system [6]. In our effective A— model the fictitious pendulum 
may start from the top (6 = 0) or the bottom (9 — tt), and come back to the same initial location, 
since its motion is limited to a single island in FigflH]) between two end points (where velocity 
vanishes). Calculation of the pulse area for this motion gives the quantized unit of 2n. 

A more physical reason for existence of solitons is that the propagating pulse in medium is not 
stable and becomes reshaped via spliting and compression. The sharpening process however does 
not last indefinitely, since the field energy density increases without bound. If the coherent region is 
maintained long enough, the sharpened pulse eventually emits PSR, the only possible process in our 
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model, thus becoming stable against PSR by formation of solitons. This picture shall be supported 
in conjunction with PSR rate computation in the next section. 
Solitons in our effective A— model are solutions of 



d v 9(r] , t) = =Fa(cos 9(r] , r) — 1) , r = t — x, rj = x , (44) 
with the boundary condition suitablle to the kink solution, 

0(77,-oo) = 0, 0(77, 00) = 2?r , (45) 

and a finite energy condition, 

POD 

dt(d t 9(x,t)) 2 < 00. (46) 



-(+) case in eq.()4*4l) corresponds to amplifier (absorber). The anti-kink, or anti-soliton, does not 
exist, because the relevant condition, 0(77, —00) = 2n , 9(r], 00) = implying d t 9 < at some time t, 
is excluded from the required positivity of field energy due to \Eq \ oc d t 9. 

The explicit form of the tipping angle 9 and the pulse strength \E 2 \ of soliton solutions is given 
in terms of a single function T(y), 

n, \ l-a 2 (x-T(t-x)) 2 

^f\E\ X ,t)\=± 2a % (48) 

4 1 + a z (x — 1 (t — x)Y 

Finiteness of energy requires a behavior of the function T(y) towards ± infinite time, 

*f|y <00 ' Sy^§ <0 °- (49) 

The function T(y) must be monotonic to give a positive definite flux. 
A simple, linear choice of T(t) oc t gives 

,9, m 2Bv(l-v) . . 

where a is defined by (I35j) . The parameter v restricted to < v < 1 here may be regarded as the 
velocity of soliton. In the v — > 1 limit, soliton becomes sharpened without bound; its width oc 1 — v 
and its peak value oc 1/(1 — v). The integrated soliton flux is 

2ttvB 4 E 2 (E 2 — E 2 I A) 

= l v ^9^9 e_gl_J_ 4 ?5 x 10 7 Wmm -2 cm x y ^ 

a 3 7 i9 
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Soliton structure: v=0.9999 
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Figure 9: Soliton sandwitch structure: field power and net emission rate, namely, 2 photon (emission 
- absorption) rate. 

(numerical value for Ba), which is independent of the target number density. 

Introduction of higher order powers in T(r), like r 3 , gives asymmetric distortion of pulses. 

Fig© shows an example of soliton structure within target. The soliton has a sandwitch structure: 
2 photon emission region in the middle is surrounded by two sides of absorbion region. This way the 
net two photon emission seen outside the target vanishes, as more fully explained in the next section 
on PSR rate calculation. These solitons are stable, keeping their pulse area of 2ir and their shape, as 
illustrated in Fig ffTOl) . The soliton size is characterized by 1/a, which is oc 1/n (inversely proportional 
to the target number density). When the soliton propagates within medium of a different, hence 
mismatched number density, its shape changes, as illustrated in Fig flTTI) . 

It should be kept in mind that solitons represent, not only the field but also, the entire coherent 
state of atoms and fields as a whole. When a soliton exits from a target end into another region 
of different environment (for example, of different target number density or of different matter in- 
cluding vacuum), a mismatch of soliton parameters occurs and the soliton becomes destabilized and 
necessarily emits PSR. This gives a simple principle of detecting a soliton and PSR at the same time. 
An obvious obstacle against this is a fast relaxation process. 

We now discuss how solitons might be dynamically created in a long target. The first step towards 
multiple soliton formation is to irradiate a strong pulse of area close to an integer n x 2tt and create 
n number of well separated pulses. Each of these pulses are candidates of soliton of area 2ir, but 
they must be reshaped. A method of reshaping would be to recast these pulses into an amplifier 
close to the ground state given by the angle 8 = ir. The reason for this is that at this angle value 
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Propagation of soliton: v=0.9999 
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Figure 10: Spatial profile invariance of propagating soliton, shown at equally spaced different times 
(the first blue one giving the original pulse). Ba number density 10 15 cm -3 and the velocity 0.9999 x 
the light velocity are assumed. 



Propagation of soliton in wrong medium: v=0.9999 
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Figure 11: Spatial profiles at different times when the same soliton as in Fig ffTUj) is put into a medium 
of a mismatched number density 10 times larger than in Figf fTU]) 
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the solution is given by eq. (l3"0"j) . thus at late times one gets 



\E 2 (x,t)\ 



x 



1 — ax5(t — x) + (ax) 2 



(53) 



where the function 5(y) is the area - 2ix, and 5(y) = 9(y) — 2n ~ {d y 6)o y. This has the same form 
as the soliton solution, eq.(H81). with T(y) = 8 /(2a) ~ y(d y 9) / (2a), if 5 is small. In Fig(jT2i) we 
illustrate a profile of nearly formed soliton, constructed by propagating a nearly 2ir pulse in a long 
medium of 9o ~ n. The pulse profile constructed this way is similar to, but a little bit distorted 
from, the soliton profile in Fig([H]). The reshaping becomes perfect when the process of PSR occurs; 
the subject of the next section. 



Nearly formed soliton: 1.04 *2Pi pulse 
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Figure 12: Spatial profile of nearly formed soliton. The net emission rate and the pulse profile after 
a long passage in the ground state amplifier are shown. #o = 0.9999 x 7r, and Ba number density 
10 16 cm~ 3 assumed. 

Two solitons infinitely separated from each other in medium propagate independently, but at 
finite distances they start to influence each other. Two soliton interaction of this kind is similar to 
the van der Waals interaction between two neutral atoms caused by induced dipoles. Interaction of 
two solitons may be calculated using the effective potential given by eq.( H3"j) in the preceding section. 
Examples of two soliton potential are shown in Fig ffTBl for various target states. The conclusion 
on two soliton interaction is that its potential has many extremal equilibrium points for all types 
of target states. As the target number density increases, the depth of potential (local) minimum 
becomes larger and its locations becomes closer to another soliton. Implications on this force nature 
shall be discussed elsewhere. 
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2 soliton potential: absorber, amplifier and dark states 




cm 



Figure 13: Effective potential between two solitons, both located in a single medium. 1st soliton 
of velocity 0.9 x the light velocity is at 5cm away from the target end, and separated from 2nd of 
velocity 0.99 x c by indicated distance. Ba number density 5 x 10 16 cm -3 is assumed, and 4 different 
target states, the amplifier (in blue), the absorber (purple), and two dark states (in amplifier and 
absorber corresponding to two signs in eqs. fll41j) and (I142p ) of 6* = 7r/2 (green and brown) are 
compared. 
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VI Theory of PSR and its relation to soliton formation 

At the outset we would like to point out that there are two kinds of PSR, the one initiated by 
a trigger and another that occurs without (or with a very weak) trigger, looking like a spontaneous 
emission at a superficial level. These may be regarded as two parts of the coherent two photon 
emission under an incident trigger pulse; the stimulated and the spontaneous parts. The stimulated 
part has a rate oc iV 7 (the number of photons in the incident pulse), and the spontaneous part oc N®. 
The rates of these two processes differ by a large factor, roughly reflecting their available number of 
photons; O[|£^o( w )| 2 / W ] (the number within the trigger pulse) for the triggered PSR and O[co 3 / (27r 2 )] 
(the number within a single photon phase space) for the trigger-less PSR. The spontaneous to the 
stimulated PSR rate ratio is numerically 

O[10~ 4 ](-) 4lQ6Wmm ~ 2 - (54) 

Although its rate is smaller, there are merits for PSR without (or with a very weak) trigger; 
easiness of two photon simultaneous detection and the presence of the PSR time delay, which can be 
used as an experimental identification method for a good PSR event. It should further be noted that 
even a very weak trigger pulse for PSR is very useful to expedite the target coherence. An optimal 
and careful choice of the trigger power is obviously required. 

We now discuss PSR with trigger of a larger rate. The backward photon emission in PSR is 
treated as perturbation to the pulse propagation. The propagation part was analyzed by Maxwell- 
Bloch equation and already solved. We are regarding the process as a whole single event, starting 
from the target triggering, coherence evolution, until PSR occurs. 

Introduce a general perturbation to the propagation problem, given by a Hamiltonian % to the 
effective A— system, which may cause a variety of transitions depending on the choice of in the 
case of PSR 1-Li = d ■ e where e is the field of a single photon emitted in the backward direction to 
the pulse. For RNPE, % ~ GpvlaVj • S e (S e the electron spin operator), with i/j the neutrino field 
of a z-th mass eigenstate. 

The interaction Hamiltonian gives rise to a perturbed amplitude of the upper level in equation 
for dc e /dt, which contains the amplitude of upper level \j), (see eq. (l88|) in Appendix I). Amplitudes 
we need to consider in this case are the unperturbed amplitude Cj(t) and the pertubed amplitude 
Scj(t), given respectively by 

c (°) ( t ) d 2^ e -^~E ]9 )t E * it^ e i(*+E je )t ECe (55) 

J U — Ejg U + Ej e 

6 c (t) d 2^_ e -^'-E ja )t e * Cg ^_ e i^'+E je )t eCe (56) 

U) — Ejg UJ + Eje 

Pulse field x perturbed backward field contain spatial functions, E oc e~ tkx and e oc e~ tk ' x , where 
(a;, k) , (a/, k') are frequency and wave number sets of respective waves. 
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A crucial observation is that at the half of two level energy difference, oscillating terms have a 
common phase, since 



-(w + E je ) = cj - E jg , (57) 

at oj = E eg /2, which much simplifies the amplitude. 

We introduce a new notation for propagating and backward emitted photon fields by 

E * e i(u,t-kx) + EQe -i(ut-kx) ^ ( 5 g) 
e * e i(u>'t-k'x) + eoe -i(u,'t-k>x) ^ ( 5 g) 

and derive the amplitude of \j) -level, cf\t) + Scj(t), by time integration, leading in the Markov 
approximation to 

cf(t) = (e^-^-^E^c^t) + e-^-^ t+lkx E d jg c 9 (t)) , (60) 

iE t 

6 Cj (t) = - C — (e^'-^-^'^ld^it) + e~ l ^'-^ t+Mx e d jg c g (t)) , (61) 

where uq ~ E eg /2 is the frequency of triggering pulse. Neglected terms are rapidly oscillating with 
x at kl ~ — k (giving two counter-propagating waves). 
This gives corresponding perturbation amplitudes, 

5c e (t) = % c g (ty^> d ^ E ^ 4 Sin{u + 2Uo)t/2 , (62) 

E c io + oj' — 2u 

Sc g (t) = ,c e (t)e-^> d ^ E ° e ° 4 Sin ^ + - 2 ^l 2 . ( 6 3) 

One may interpret eo and as annihilation and creation operators of backward emitted photons 
according to quantum field theory. 

Equation f )63|) describes paired emission of two photons, counting the stimulated part of amplitude 
(oc WNj) of incident pulse, and ( |62l) describes an associated process of paired annihilation of two 
photons. The reason why the paired annihilation occurs is that in medium under propagating field 
some atoms may be in the ground state surrounded by ambient two modes of photons. This amplitude 
is Sc g , the perturbed ground state amplitude. 

When one computes the probability of PSR, one first sums amplitudes including e ±l ( k+k ') x over all 
atoms in a coherent medium. For this purpose we introduce local atomic amplitudes c;(x, t) ,i = e,g 
and divide the entire coherent region < x < L into cells of size Aa;. We require 1/k <C Ax <C L 
for the cell size. Under SVEA, namely the assumption of slow variation of local amplitudes over the 
wavelength scale, the amplitude summation within a cell gives 



X+AX dy Ci (y, t)e^ k+k '^ ^ Cl (x + ^, t)e ^ + y)A x /2 ^Mk + k')/\x/2 ^ 
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The last factor, when squared, gives the Dirac delta function in the form, Ax2^S(k + k'), in the large 
Ax limit (of Ax 2> the wavelength). The factor Ax in front is cancelled by the amplitude squared 
|e | 2 of a single backward photon. 

The factor 27r5(k + k') implies the momentum conservation working, independently of the actual 
finite, but a large Ax. Our basic ansatz in the spirit of SVEA is to take for Ax its largest possi- 
ble value, the target size L. Furthermore, we replace for ease of numerical computations the fast 
oscillating function by a more smoothly varying, yet globally correct, function; 

sin(fc + fc')L/2 \ 2 1 , * 

k + k' J {k + k') 2 + 4/{L 2 ) ' 1 ' 

We may, somewhat arbitrarily but realistically, change in rate estimates the factor L here by its frac- 
tion, which gives larger rates than given here. We may thus regard our procedure of the replacement 
Ax — > L as an underestimate of more realistic rates. How much rates are underestimated actually, 
however, must be verified by a more laborious numerical simulation, assuming discretized sites of 
atoms. 

Rate of the net emission oc R s has two types of contributions; emission for i? 3 > from atoms 
populated more in the excited state, and absorption for R 3 < from atoms populated more in the 
ground state. The absorption that occurs within target medium cannot be experimentally measured, 
and one measures the positive emission rate at both target ends. (If the effective rate becomes 
negative, no emission is measured since it means an inward emission into the target inside.) The 
spacetime dependence of net rates is thus given by 

\E, *|» nRM t) J(x, «)(4 l , + ,,_ 2 ; o W > 2 (^f) 2 ( ¥WTW ) • (66) 

where the pulse related factor I(x, t) is to be given below in eq. fl73|) . We may define the net spectral 
rate (net probability per unit time which can also become negative) as 

, o, s / wWoni /sin (a; — u')L/2\ 2 , , 

\El\WunRz{x,t)I{x,t){^-^) 2 *A (—± U-) , (67) 

Cj c \ U) — it) J 

by the short time average using the well known formula, 

.. 4 /sin((u; + u/-2u;o)i/2)y n , , , aoS 

hm - — ^ ' ' ' = 2tt5(uj + J - 2lo ) . (68) 

t^oo t \ to + to' — 2ujq J 

The large time limit is valid since oscillation in time is very fast compared to polarization development 
under all practical situations we consider. 

Convolution with the frequency distribution of the input pulse is now necessary. We assume a 
Gaussian power spectrum of the form, 

771/ r (^-^o) 2 ! , Rn , 

F(uj;ujo,6) = — -^-exp[ ], (69) 
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with the frequency width S(— 0(50 ~ l)GHz for commercially available laser). The wave number 
integration, along with the convolution, gives 

J™ du jH dk' j^^^ 5{u + \k'\ - E eg ) ~ ^F(E eg - J ; c , 5) , (70) 

leading to the differential spectrum d 2 T /dui'dt (having the dimension of rate , 1/time) of the backward 
photon of energy u/; 

a 2 2t/ /o ,2 7je7j 3 ^e 3 7T E eg - u' e 2 (ajg) exp [- (gl -y"*) 2 ] exp [_ ggg] 

87r n v(37r) ^^;~ (^-^ + 4/^ 7^ 7^ rs(x ' * ; u } J(x ' * ; w } ' 

(71) 

1 cosh(o;x sin ) cos #o(l ~ cos $o cos $) ± sinh(o;x sin 9 ) cos # sm $o srn $ — ( cos $o — cos $) 



y/l + 7 2 cosh(aa;sin #o)(l _ c °s #o cos#) ± sinh(o; sin #0^) sin^o sin 9 — cos 9 (cos 9 — cos 0) 

(72) 

/(:M;a/) = = — = (73) 

cosh(o;a; sin $o)(l — cos#o cos 6?) ± sinh(aa; sin 9 ) sm6 sin 9 — cos 9 (cos $o ~~ cos 6?) 

For the input pulse we assumed a time structure of the Gaussian form, characterized by a width A 
and centered at t = x. Two cases of =F correspond to the state of targets; absorber (amplifier). 

Note that we used neither the conventional Gaussian width given by the variance, nor the half 
width at half its maximum. Our width is \pl x the Gaussian and 1 / Vhi 2 x the half width. 

For the Ba D-state, the basic unit of differential spectrum for the target number density n = 
10 16 cm~ 3 and the pulse power e 2 , = 10 6 W muT 2 is (under the h — 1 , c = 1 unit) 

87r 2 (37r) 2 (10 16 cm- 3 ) 2 cm 2 GHzMHzl0 6 Wmm~ 2 = 5.7 x 10 19 HzeV, (74) 

to be multiplied by 

( n V E eg — oj' e 2 e^[- gfggg ]jy 

4o 16 cm- 3j cm 3 (2a/ - E eg ) 2 + 4/L 2 10 6 Wmm- 2 ^5 

exD r_(£z£)f.i fa 

xr 3 (u' , x, t)I(u>' , x, t) P[ A2 J , (75) 

a/7TA 

where all energies should be given in the eV unit. Two cases of +(-) correspond to the absorber 
(amplifier). 

The time dependence of rate d 2 T/du'dt is due to the rapidly varying pulse shape. A more practical 
measure of event rate is time integrated event number per a shot of pulse; 



"** r fr' - (76) 

du'dt y 1 

From the measurement point of view, only two locations are of practical importance; x = the 
target left end for measurement of the backward emitted photon of PSR, and x = L the right end 
for measurement of the forward emitted photon. 
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The frequency spectrum of backward photon is sharply peaked at the middle point of level 
spacing, uj' = E eg /2, if the trigger frequency is tuned to this value, u = E eg /2. A quantity of 
practical importance for the overall rate is the frequency integrated (time dependent) rate given by 



o 



d*iy ; s,t) _ r eM- {u '- E t W0) W E eg -u> . 

dU dJdt - I dU (2u>-E eg y + 4/L* F{UJ (7?) 



where F(u' ; x, t) is slowly varying in frequency u/. We may approximate this frequency integral by 
taking correct behavior in the two limiting regions, 5L ^> 1 and 5L <C 1, and smoothly interpolating 
in the interval. At the tuned point Uq = E eg /2, one may adopt the following approximation, 

v/i<S (2u,'-B eJ )2 + 4/L2 1 ' ' ' ^+2LS ' 1 ' 

We checked that this approximation is valid to an accuracy of ~ 15% level. 

Assuming that the pulse shape variation is slow, one may introduce a slowly varying, time de- 
pendent rate by taking the short time-average and frequency-integration, to obtain a total net event 
number (emission - absorption event) per a shot of pulse as 

Af(x, t) = 3.2 x i "(— J . 2 —,rA x, t)A x, t)^^- . (79) 

10 lb cm 6 lCrWmm cur 2 2 y'n + 2Lo 

The last factor for the pulse duration of 5 = 45GHz and L = 3 cm is ~ 0.5. 

The more effective and useful event rate (number of events per unit time) in actual experiments is 
given by this event number per a shot divided by the repetition cycle time r r of excitation and trigger. 
For instance, if the pulse repetition cycle r r is 1 msec, a practical experimental rate is M /r r = 10 3 A/" 
Hz. We caution that all event numbers shown in Fig( fl4|) ~ Fig( l23|) are events per a shot of pulse, 
and one has to multiply l/r r (depending on experimental setup) for effective and more realistic rates 
per unit time. 

We now exhibit several figures to illustrate physics of analytic results given here. These PSR rates 
are computed, using the rate formula and analytic solutions of the pulse propagation for r 3 , I^. The 
first figure Fig (fT4j) shows the initial angle 9q dependence of the backward rate (measured at x — 0) 
and the forward rate (at x = L). Except at soliton formation discussed later, the forward-backward 
asymmetry in rate is not large (typically < 15%). 

An example of Ba backward photon spectrum is shown in Fig( !T5|) . In this computation we took 
the Gaussian frequency distribution of input laser, its width given by 45 GHz. Both spectral shape 
and rate are indistinguishable up to the passage time of pulse at the target end, but the rate rapidly 
decreases much beyond the passage time. It thus becomes important to devise a fast recycling scheme 
for excitation and trigger of the target, in order not to wait for a null result. 

The laser power dependence of rate is shown in Figf JT6|) . Fig (fTT|) shows dependence on the target 
number density dependence, indicating that PSR may be detectable for Ba number densities as low 
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events Ba fb rate vs angle 




Figure 14: Angle (#o) dependence of forward (blue) and backward (purple) rates measured at the 
passage time of the target end (taken 3cm). Power 10 6 W mm -2 , number density 10 16 cm -3 , laser 
duration 3 ns assumed. 
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Figure 15: Spectrum rate at 3 different times, (0.1, 1, 5) x the passing time of the target end. Power 
10 6 W mm -2 , number density 10 16 cm -3 , target length = 3 cm, laser duration 3 ns, 9 = tt /A assumed. 
Rates monotonically decrease in this range of time sequence. 
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Figure 16: Rate at the arrival time of the target end as a function of laser power. Target number 
denisty 10 16 (in blue), 10 14 cm 3 (in purple), target length = 3 cm, laser duration 3 ns of 45 GHz, 
and 6q = tt/4 are assumed. 



as 10 4 cm~ 3 or even less, by taking into account a fast repetition cycle time of, for instance 0[/xsec] 
(O[10 6 ] times the event number of Fig(ll7|) for the rate per second), whose precise value is determined 
by experimental conditions and not by a theoretical calculation. Dependence on the number density 
is roughly oc n 2 . This may open another interesting possibility of detecting PSR in alkhali earth ions 
such as Ca + , Sr + and Ba + , which have the A— type level structure. The actual PSR rates, calculated 
by using experimentally known level spacings and decay rates for these ions, are somewhat smaller, 
as illustrated in Fig (fT8l) for Ca + ion. Maintaining coherence in trapped ions might however be easier 
than in the gaseous phase. (Even production of the crystalized ion has been achieved.) A fast 
repetition cycle of excitation and trigger laser of order 1 msec might be quite sufficient to obtain 
detectable rates under good envirornment of the ion trap. 

One might have a suspicion that PSR is not detectable, because a single photon SR rate is 
always larger than the two photon PSR rate due to a larger spontaneous emission rate, when both 
considered as elementary processes. This is not always true, because what should be compared in 
triggered experiments is the triggered PSR time vs the delay time of trigger-less SR, and these two 
times scale with the target number A^ as oc N~ 2 vs oc N~ l . With a sufficiently large N, the triggered 
PSR occurs before the delayed SR occurs. 

Incidentally, (J=0 — > J=0 ) PSR transitions, certainly present for alkhali earth atoms (Sr, Ca etc.) 
and Yb, are very interesting from the point of realizing excellent quantum entanglement, because two 
emitted PSR photons are well entangled in their angular momenta, back to back emitted directions, 
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and identical photon energy, all to good accuracy. 
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Figure 17: Rate at the arrival time of the target end as a function of Ba number density. Two power 
values 10 6 (in blue), 10 2 W mm -2 (in purple), target length = 3 cm, laser duration 3 ns, 9 = it /A 
assumed. 

Fig(Tl9l shows a spectrum rate at a mistuned frequency of the input laser. Note that for an 
infinitely long medium the momentum and the energy conservation forces PSR spectrum to have a 
5— function like peak at the half energy. In practice, the target has a finite length and violation of 
the momentum conservation leads to a small tail away from the peak location of u = E eg /2. The 
amount of suppressed tail contribution depends much on the frequency distribution of the tail part of 
irradiated laser. A large suppression factor seen here is due to the Gaussian frequency distribution, 
in this case of width 45 GHz. This large suppression is encouraging from the point of enhancing the 
signal to the background ratio of RNPE/PSR. We shall have much more to say on this in the last 
section when we discuss prospects for RNPE. 

We now demonstrate that solitons are stable against two photon emission. As discussed above, two 
photon emission associated with pulse propagation accompany simultaneously two photon absorption, 
since it may be induced by surrounding field. Thus, the net emission rate is in proportion to -R3/, 
the product of population difference and the power of propagating pulse. For soliton solutions the 
time integrated rate at the target end is 




(80) 
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Figure 18: Rate for Ca + ion vs number density. Power 10 6 W mm 2 , target length = 3 cm, laser 
duration 30 ns of 1 GHz width, 9 = it/ 4 assumed. 
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Figure 19: Photon sepctrum in blue when the input laser frequency is mistuned by 10~ 4 away from 
the middle energy of E eg /2. For comparison the original spectrum at the tuned frequency is also 
shown in purple. 
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From Fig([H]) one sees that the emission region in the central part of the target is sandwitched by 
two absorption regions, which gives a balanced vanishing net rate. This result holds at any target 
point. In more general target states excluding solitons, the integral of this product oc cos8d t 6 gives 
a difference, sm0(t = oo) — sin#(t = — oo), which is non-vanishing. 

An ideal method of observing PSR might be creation of many solitons at the first stage, and their 
artificial destruction by controlled means for detection of PSR photons at the second stage. 

The forward-backward asymmetry is expected to be large, both immediately prior to soliton 
formation and immediately after their destruction. We plot in Fig( )20|) time evolution of the asym- 
metric rates for ~ 27r pulse. This asymmetry may be used to detect PSR itself under large symmetric 
backgrounds. 



FB asymmetric rates: 2Pi pulse 
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Figure 20: Forward-backward asymmetric PSR rates. The backward rate in blue is the measured rate 
at x = 0, and the forward rate in purple at x = L (the target length 10 cm taken here). The rate here 
is (emission - absorption), hence the negative value region gives zero measurement. 6q = 0.9999 x n, 
and Ba number density 10 16 cm~ 3 are assumed. 



VII Outlook for radiative neutrino pair emission (RNPE) 

We shall briefly sketch prospects towards our goal of the precision neutrino mass spectroscopy, by 
providing rate, spectrum and S/N. Details of rate computations, including effects of all six thresholds, 
mixing angles, Majorana vs Dirac distinction [I], and Majorana CP phases [5] shall be presented in 
a separate publication. 

Atoms ideal for PSR may not be appropriate for RNPE detection, due to two reasons: (1) the 
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RNPE process |e) — > \g) +7 + 1/^ , i, j = 1, 2, 3, requires a large spin flip amplitude [I] (similar to Ml 
transition for the electronic part of transition matrix element) between the intermediate state and 
either of the initial or the ground states, while a large PSR rate may require larger El transition, 
(2) PSR photons might become a serious background against RNPE, which means that smaller 
PSR rates are better for RNPE detection. To the best of our knowledge, Xe atom in solid form 
is among the best candidate atoms for RNPE. The candidate metastable state is the first excited 
state, 5p 5 ( 2 P 3 / 2 )6s 2 [3/2] 2 , a J = 2 pair state of electron and hole [5]. Other candidates might be 
metastable states of spin configuration (of two electrons) different from the ground state, commonly 
seen in diatomic molecules such as O2, an object worthy of serious consideration. 

Experiments can be performed using typically more than three lasers of different frequencies; more 
than two for excitation to the metastable state and another for the trigger of RNPE. The trigger 
frequency u (different from E eg /2, hence mistuned for PSR) is reset each time for measurements 
at different photon energies of the RNPE continuous spectrum of Fig fl2Tj) . Hence it is desirable to 
use frequency tunable lasers for the trigger. This way there is no Gaussian tail suppression at each 
detected spectral point of RNPE photon, while the background PSR is suppressed by the Gaussian 
tail factor due to a mismatch away from u> = E eg /2. Moreover, the energy resolution of RNPE 
photon is essentially determined by the precision of triggering laser frequency, and not by detected 
photon energy resolution. This is a key for success of the precision neutrino mass spectroscopy, which 
must resolve photon energies at the /xeV level. 

Theoretical estimates readily give coherent RNPE rates as large as, of order 0.1 events per pulse for 
Xe atoms of number density of order 10 18 cm -3 , which may be realized in solid matrix environments. 
With a repetition cycle of 0[1] msec interval, this gives a detectable effective rate of O[100] Hz. 
We show the continuous single photon energy spectrum of Xe macro- coherent RNPE in Fig lpTj) . 
The sharp rise at the threshold is characteristic of the three-body decay under the momentum 
conservation, as is familiar in the \i decay. The increasing rate towards the low energy photon, given 
in the blue curve of Fig fl2Tl) . is due to the flat frequency dependence of the peak intensity of the 
trigger laser, and the decreasing rate in purple is due to oc u 4 peak intensity dependence more akin 
to the phase space of the 3-body spontaneous decay of elementary particles. The low energy side 
of the photon energy spectrum is thus sensitive to the frequency dependence of the trigger laser 
intensity, and is inevitably tied to experimental apparatus used. The rate scales as n 2 V, with the 
target number density n and the volume V of coherent region, as shown in Fig (l2~2l) . 

The problem against the precision neutrino mass spectroscopy is not the rate itself, if a sufficient 
number of target atoms is prepared. Rather, the serious problem appeared to be in the signal to 
the background ratio (S/N) where the main background source is the physical process of two photon 
emission, in particular N 2 enhanced PSR. We discuss this problem shortly. Higher order QED 
processes that could sneak into our photon energy region might appear problematic, but they are 
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Compared RNPE rate per pulse: n=10 A 18 cm A -3, l=10cm 

events 




photon energy in eV 



Figure 21: RNPE spectrum rate (event number per pulse) starting from the neutrino threshold 
of the pair of mass = 50 meV. Excited Xe number density = 10 18 cm -3 . Two different frequency 
dependences of the laser peak intensity, the flat one (in blue) of magnitude 10 6 W mm -2 and the 
one oc u A (in purple) are compared. Time duration = 3 ns. Complications due to mixing and phase 
factors, and to other thresholds are all ignored. 




Figure 22: RNPE rate (event number per pulse) vs the target number density of pair emission of 
m= 1 meV evaluated at m=50 meV pair threshold. Laser intensity 10 6 W mm -2 and time duration 
3 ns are assumed. 
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actually negligible if they have (spontaneous) photon emission rates smaller than the elementary 
decay rate of metastable (lifetime > 0[1] msec in our standard) state, since in actual experiments 
a cycle of measurement for RNPE is terminated much earlier than the lifetime of metastable atom, 
and the target preparation is recycled. The single photon SR, entirely outside the photon energy 
region of our interest, is negligible also in rate by the choice of trigger frequency for RNPE. 



Xe spectrum: RNPE and PSR, 10 A 17 cm A -3, 10cm, 1 GHz 

events 




4.14 4.14 4.14 4.14 

Figure 23: Rate of mass 1 meV pair emission (in blue) and the background PSR rate (in purple) very 
near the neutrino mass threshold. Excited Xe number density 10 17 cm -3 of target length 10 cm and 
laser intensity 10 6 W mm" 2 of frequency width 1 GHz of time duration 30 ns are assumed. When the 
blue curve is above the purple curve, S/N > 1. At the farthest left is the threshold of the pair mass, 
(8 + 1) meV, and at the farthest right is that of the pair mass, (1 + 1) meV, the entire abscissa 
distance in this figure being separated by a very small energy interval ~ 4.6/ieV or ~7.1 GHz, hence 
their precise values are not shown. 

As shown in the present paper, the S/N problem caused by PSR is much relaxed and the value of 
S/N is enormously enhanced by two effects; (1) mistuned trigger frequencies for PSR and (2) soliton 
formation. We first discuss the first issue. S/N is increased by the choice of trigger frequency set 
at uj < E eg /2 — 2m\/ E eg of the lightest neutrino pair emission for RNPE (we are ordering neutrino 
masses according to > rri2 > mi). There is a mismatch of energy from the largest rate point at 
u = E eg /2 of PSR by the amount 2m\/ E eg . If one may assume the Gaussian frequency distribution 
for the trigger laser, the PSR background rate is suppressed by a large factor of exp[— 4mf/(E^ g S 2 )}, 
with 5 the Gaussian frequency width. A numerical factor of practical importance is ^jE eg 5/2 ~ 
1.2meV-y/ 8/ GHz. It is thus important to have the frequency width of sub-GHz in order to explore 
the sub- me V neutrino mass range. 
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In Fig (l2"3"l) the spectral RNPE rate of assumed neutrino mass threshold is plotted, along with the 
background PSR rate, taking for the input laser 1 GHz frequency width of Gaussian distribution. It 
is estimated from this and similar figures of the pair mass range, mi = (0.5 ~ 1) ,m 2 = (8 ~ 10) 
meV, that this quality of laser can explore the neutrino mass of order 1 meV. With an even narrower 
width laser, one may hope to explore the sub-meV region and detection of the relic neutrino [9]. 
This, however, rests with the Gaussian nature of the tail part of the laser frequency distribution, 
which must be carefully examined from experimental points. 

Soliton formation gives another great merit in which the background PSR is blocked, completely 
in ideal situations, and may open a new path towards a controlled RNPE experiment. In practice, 
only a partial blocking can be expected due to various kinds of relaxation processes. Naively, one 
would expect that the ratio of the relaxation to the soliton formation rate is a relevant quantity 
for the partial blocking. The soliton formation rate is however difficult to estimate without detailed 
experimental design and experimental R&D works. Study of decoherence in solid environments by 
means of detailed numerical simulation and R&D works is also very important for realization of 
the neutrino mass spectroscopy. 
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VIII Appendix I Derivation of effective two level model 

• Atomic system 

The state vector of an atom can be expanded in terms of the wave function, 

m)) = E c ^) e ~^'> + c ^y~ lE9t \9) + Ce(t)e-^\e) . (81) 
j 

c a (t) are probability amplitudes in an interaction picture. 

The atomic system may interact with fields. The electric field e(x, t) that appears in the 
Hamiltonian via El or Ml transition is assumed to have one component alone, namely we ignore 
effects of field polarization. This is a valid approach under a number of circumstances. One 
then decomposes the real field variable e(x, t) into Fourier series, e %UJt times a complex envelope 
amplitude E (x,t), and its conjugate, where E (x,t)e tkx {k = wave number) is assumed slowly 
varying in time, 

e(x, t) = E (x, t)e iwt + E*(x, t)e~ iujt . (82) 
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The Schrodinger equation for a single atom 



t^(t)) = (H + de)\m), (83) 



gives the upper level amplitude Cj(t). Using 


''dt 

one has 



£(j\m) = (mo+de)m)), m 



rip. 

l lt 6 ~ lEjt = ( d ^ e ~ lEet + d ^ 9 e- lEgt )e . (85) 
This can formally be integrated to 

c 3 (t) = ~ij*dt (dj^x^c^t'y^-^ + djgeixrfcgWe'W-W 

= -i J dt' (d je c e {t')e iE ^' + d jg c g (t')e iE ^ (E G (x,t')e iut ' + E*(x,t')e~ iult '^ , (86) 

(the initial condition Cj(0) = is assumed). 

We extend the system of a single atom to a collection of atoms, regarding variables c a (x,t) 
as functions of two variables, x and t. The population and the coherence of the entire atomic 
system is described by the density matrix elements, R a p(x, t) , a(3 = e, g,j, which is the squared 
product of wave functions (bilinears in c e ,c g , Cj and their conjugates) and the target number 
density n(x). The macro- variables are defined by 

Rap{x, t) = n(x)c* a (x, t)cp(x, t) . (87) 

For a notational simplicity we often omit space coordinate (x) dependence of probability am- 
plitudes and write these simply as c a (t) instead of more proper c a (x,t). 

Markovian approximation and SVEA 

The basic strategy of deriving equations for the lower two level amplitudes c e , c g in a closed form 
is to eliminate atomic variables related to the upper level amplitude Cj. This is essentially done 
by neglecting a long-time memory effect (the Markovian approximation) and making slowly 
varying envelope approximation (SVEA) in the terminology of [3]. The idea is to replace 
dynamical variables, c e (t'),c g (t'), E(x,t') in the integrand of eq. (l86|) . by their values at time t, 
neglecting all the past memory effects. This gives 

/ 1 — e ^+E 3g )t i _ e -i(u-Ej B )t \ 

+d " C ' { U+Ei , B » ^E— E "°) • 
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which is inserted into equations for the lower levels 



^l = ^J2d ej e(x,t) Cj (t)e iE ^, (89) 



dt 



Unitarity, namely the probability conservation, given by 

!(|c e r + |c,| 2 + 5>| 2 )=0, (91) 

j 

does not hold when limited to two amplitudes c e , c g alone, namely in the Markovian approxi- 
mation. However there exists an effective conservation law that holds for the two level system, 
given later by (11141) . There seems some misunderstanding on this point in j3]. 

We further neglect rapidly oscillating terms assuming the nearly, but not necessarily exactly, 
tuned condition, u ~ E eg /2, which amounts to 

= i (^ ee \E 2 \c e (t) + » eg (E* ) 2 c g {t)e-^- E ^) , (92) 

= * (»Mc 9 (t) + li ge Elc e (t)e^-^)t) , (93) 

d 2 E ■ cP E '■ 

/^ = 2 E EI-oj^ ^^^/-I 2 ' (94) 

j je j jg 

,, — djedjg _ dj e dj g , . 

j j 

where E c = |(-% s + Ej e ) , 5u = u — ^ as in the text. All of dj a , \i a p are taken as real. 

A mistake in [3] is that \x eg = /i ge is assumed even at 5u ^ 0. This equality holds only at 
the middle point of frequency u> = E eg /2. Thus, results of [3] away from this tuned frequency 
should be taken with skepticism. 

Writing this equation in a matrix form, 



d 
dl 



-iU ^ , (96) 



C 9 I \ C 9 



, Vee\E*\ e-«*«- E ^% g (E* o y 
1 1 e^~ E ^^ ge El ^ gg \El\ 1 ' 



one finds that the effective Hamiltonian % becomes hermitian only by neglecting 5u/E c terms 

m fi e g ; l^ge- 
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We introduce symmetric and anti-symmetric functions of Sou = u> — E eg /2 as 

, \- d je d jg E c \- d je d jg 5u 

" = 1^E[^> " =l^Ef^> (98) 

3 3 C 

Heg = /i + + /i" , Hge = fi + ~ fi~ , (99) 

and write 

% = %+ + %-, (100) 

" ' e*(a-^V£g V 99 \E 2 \ ) ' 

?T = , , . 102 

I _ e i(2 W -E eg )t At - jE ;2 o 7 

"H + is hermitian, while "H~ is anti-hermitian. The anti-symmetric piece %~ vanishes at the 
middle point of 5u = 0, since H~ oc pT oc <5w, hence is small except where 5u is large of 0[i? c ]. 
We shall ignore effect of H~ and assume \i ge = /i eg in the main text of the present work. 



• Generalized Bloch vector and its dynamical equation 

The 4 component Bloch vector is defined by 

R = n(\c e \ 2 + \c g \ 2 ) , (103) 

R 1 = in{c g c e ^ - c\c g e~*) , (104) 

R 2 = - n {c g c e e^ + c* e c g e- 111 ) , (105) 

R 3 = n(\c e \ 2 - \c g \ 2 ) , (106) 

t] = (2uj- E eg )t - 2kx + 2(p , (107) 

where we assume a standard form of field, 

e(x, t) = e (x, t) cos(ut — kx + ip) , (108) 

with k = ±u. The real amplitude €q(x, t) and the phase ip(x, t) are assumed both slowly varying 

in time and in space. Note the relation of real and complex field, = 4|E | 2 (E is complex 
including e 1 ^*^). 
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The generalized Bloch equation is given by 

d_ 
dt 



^Ri = (^j^e 2 + (2co - E eg + 2^)) R 2 + ^elR 3 - ^e 2 R , (109) 



d_ 
dt 



R2 = ~ (f^^lel + (2c - E eg + 2^)) R 1 , 



(110) 



m R9 = ~ t 2^ Ru 

^Ro = ^lRu (112) 

,,+ _ \^ d e jd gj E c (t\<3\ 

* " \ m - (- - Ej2Y ' (113) 

proved by using the Schrodinger equation for c e ,c g . Note that in RHS of these equations there 
is no phase factor like e ^t-2ut+E eg t) _ 

The conservation law is extended to the 4- vector; 

^(Rl + Rl + Rl + R 2 3 )=0. (114) 



Only at the middle point of u = E eg /2 this conservation reduces to the usual type of conser- 
vation for the 3- vector norm, §[(R 2 + R 2 + R\) = (the assumption taken in [3]), since in this 
case n~ = and one has separately d t R% = 0. 

Polarization of medium 

Polarization vector is defined by 

P = n(x)(ij(t)\d\^(t)} = n^2(d ej c* eCj e- iEjet + d gj c* g c j e- 4E *> t + c.c) , (115) 

3 

for which we eliminate Cj using eq. (l88|) . 

The polarization can be decomposed into the in-phase cos(ut — kx + if) and the out-phase 
sin(u;£ — kx + ip) parts; 

P = n(x)(fi ee \c e \ 2 + fi gg \c g \ 2 - {ii eg c g c* e e~ ir> + fi ge c*c e e lv ))e (x,t) cos(ut - kx + ip) 

-i(fi eg c g c* e e~ tv - iig e c*c e e iv ))e (x, t) sin(u;t - kx + <p) 

_ ^ee + ^99 n _|_ /^££ — ^-R 3 + fj, + R 2 )e (x, t) cos(ut — kx + (p) + fx + Rie (x, t) sin(u;t — kx + cp) . 
2 2 

(116) 

0[/i~] terms do not contribute to hermitian polarization P. 
The Maxwell equation 

, d 2 d 2 , d 2 P 



' dx 2 dt 2 dt 2 
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;ii7) 



gives for envelope amplitude and phase variation under SVEA 

{d t + d x )el = u^elR x , (118) 
ft + d x )2ip = u (-up+B* + ^^IR 3 + t^EMln^ . (n 9) 

In order to simplify equations, we introduce new variables by 

, _ R 2 - jR 3 _ 7 i? 2 + R 3 

2 ~ /t~t — 2 ' 3 ~ /t~t — 2 ' ^ ^ 

u R = ^±^^el (121) 
7 = ^Zi^, Q = 2u-E e9 + 2^. (122) 

The basic set of equations is given by 

d t Rt = u R R' 3 + —JL=(R' 2 + 7 R' 3 ) - V-e 2 Q Ro , (123) 
v/1 + 7 2 2 

^ = —r2=Ri , (124) 
V 1 + 7 2 

™o 

d t R' 3 = -u R R 1 - J. - Ri , (125) 

1 + T 



from which it follows 



d t R = ^-e 2 ^ , (126) 
(d t + d x )u R + kur = uj{i + urRi , (127) 
(d t + d x )tt = 2uy + R l VL , (128) 



dt {R 2 + Rj + {R' 2 f + (R' 3 ) 2 ) = . (129) 



We have not included atomic relaxation effects given by parameters, T\ ,T 2 ,T 2 * pQ. 

One can consistently take Q = 0, namely d t (p = E eg /2 — oj. This adjustment of field phase is 
assumed in the text of this paper. 

For both yT = and Q = the system further simplifies to eqs. (jl5l) . (1161) and (jl8p in the text 
by an appropriate choice of variables. 



IX Appendix II Details towards construction of analytic solutions 

We start from discussions that lead to introduction of the tipping angle 9(x,t), eq. ffTT|) . related 
to the Bloch vector component by R3 cx cos 8. Another important relation f fl8|) . dtO oc \E 2 \, gives a 
physical content of the area function 9, relating it to an integral of the pulse flux \Eq\. 
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Our method for solving non-linear partial differential equations of two independent variables (x 
and t) is to integrate in one variable t and replace integration constants obtained this way by functions 
including another variable x. The method works for our problem of one mode propagation, but it is 
not a general mathematical method. 

At finite time t we allow the integration constant t ,6 of Od solution (solution without space 
dependence), eq. (l27]) . to vary in spacetime according to 

to -> T{t -x), O ^ 6 p {x) , (130) 
{d t + d x )9 + a{cos9-cos9 p ) = 0, (131) 

noting a trivial equality, (d t + d x )T(t — x) = 0. Hence, solutions are written in terms of two functions 
to be determined by the initial and the boundary data, 

„, x cos 9Jx) cosh (a(t — Tit — x) ) sin 9 V ) — 1 /„„„x 

6(x , t) = arccos p ; \ y \ [ } *-L— , (132) 

cosh [oc[t — 1 (t — x) ) sin 9 P ) — cos 9 p [x) 

Ur U t s = g t9 = -a{l-d t T{t-x))sm 2 e p {x) ^ 
cosh (a(t — T(t — x) ) sin 9 P ) — cos 9 p (x) ' 

where 9 p (x) ,T(y — t — x) are yet to be determined. 

The following, somewhat complicated steps leading to eqs. (I137P and (11401) are processes of how 
the initial and the boundary conditions determine the unknown functions, 6 p (x) ,T(y). 

The given boundary data at some spatial point x — 0, the target end at which laser irradiation 
takes place, and the initial data at t — are 

6(x = 0,t) = 6(t) , 6{x, t = 0) . (134) 

These are related to the variable un(x,t) by using d t 9(x,t) = un(x,t), 

0(0, t) = ! dt'u R {0,t') , (135) 

J — oo 

0(x,O) = [ dt'u) R (x,t') . (136) 



Note that two data (11341) are independent. 

We may solve for the unknown function T(r) using the boundary condition, to get 

r^f x 1 , 1 — COS 8n cos 0{t) „ , . . . 

TOr) = r -arccosh , 9 = 6JQ) . 137 

V ; a sin O cos0 o -cos0(r) pK J K J 

The solution for 6 p (x) is obtained from the initial condition, given by 

<9f0(— x)(cos 9 P — cos0(x, 0) ) = (cos 0o ~ cos 9(—x) )ur(x, 0) . (138) 
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This equation together with eq. (1133j) calculated at t = 0,x = and (jl37j) . gives, with p, 



(flgg - [lee) 2 + 4// 2 e /4, 

cos 0(0,0) = ±cos#(0), 0(0)=/}/ dy\e 2 Q (y)\, (139) 

J — oo 

+ (-) corresponding to amplifier (absorber). The second equation of fl 1 3 9 [) is derived from the area- 
intensity relation, d t 6 oc | e§ | , for the incident pulse. 

If the major part of input laser is still far away from the target end of x = L, 6(0) m 0, and 
6(0, 0) ~ (it). The amplifier case of 6(0, 0) ~ corresponds to a physical situation in which medium 
is excited by other lasers, while the absorber case of 6(0, 0) ~ tt to medium in the ground state. 

Equation (1 1 3 3 p calculated at other points of t = gives 

cos6(x,0) = cos W - ^(^(^l^)^^^^ • (140) 

Since T(—x) here is already given in terms of 6(—x) by eq. (I137j) . this equation determines co$6 p (x) 
in terms of the initial data cos6(x,0). We find it possible to construct solutions of these equations 
only when 6 p (x) = 6q (spatially constant). 

Explicit form of solution is then (0 < 6q < tt) 

sin 2 6 \el(t - x)\ 



\Ei(x,t)\ 



cosh(ax sin 6 )(1 — cos6 cos 6) ± sinh(aa;sin 6q) sin# sin 6 — cos #o(cos 6q — cos 6) ' 

(141) 

cosh(aa;sin^o) cos^o(l — cos#o cos 6) ± sinh(ax sin 6q) cos6q sin 6q sin 6 — (cos#o — cos#) 
cosh(ax sin6 ) )(l — cos 6 cos 6) ± sinh(axsin 8 ) smd sin 6 — cos ^(cos^o — cos 6) 

(142) 

ft 

R 3 (x,t)= cos 6(x,t), (143) 

a/1 + 7 2 

with variable dependence given by x explicitly and t — x in 6(t — x). These solutions are given in 
terms of the strength of input pulse, 

f t— X 



6(t-x)=pl dy\el(y)\. (144) 

oo 



The constraint cos 2 6(x,t) < 1 is satisfied for any #o- 

Note that 6 is monotonically increasing function of its argument. The solution 6(x,t) is valid 
only for d t 6 > for the positivity of the pulse strength. 

We do not know whether these solutions have complete generality and no other solutions exist, 
but they seem to be adequately general for our purposes. The initial data 6(x, 0) at x > consistent 
with these solutions is given in terms of the input area at negative arguments 6(—x), and they are 
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independent of the boundary data 8(0, t) at t > given by the same area of positive arguments. One 
might consider the situation of 8(—x) = at x > 0, which gives 

cos 6( X ,0) = ^^ocosh(axsing ) + l 
cosh (ax sin o^) + cos #0 
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